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Abstract 

An analytical and experimental Investigation 
of the effects of mlstunlng on propfan subsonic 
flutter has been performed. The analytical model 
Is based on the normal modes of a rotating compos- 
ite blade and a three-dimensional subsonic unsteady 
lifting surface aerodynamic theory. Theoretical 
and experimental results are compared for selected 
cases at different blade pitch angles, rotational 
speeds and freestream Mach numbers. The comparison 
shows a reasonably good agreement between theory 
and experiment. Both theory and experiment showed 
that combined mode shape, frequency, and aerody- 
namic mlstunlng can have a beneficial or adverse 
effect on blade damping depending on Mach number. 
Additional parametric results showed that alternate 
blade frequency mlstunlng does not have enough 
potential for It to be used as a passive flutter 
control In propfans similar to the one studied 
herein. It can be Inferred from the results that 
a laminated composite propfan blade can be tailored 
to optimize Its flutter speed by selecting the 
proper ply angles. 

Nomenclature 

[A] matrix of Influence coefficients, 

Eq. (5) 

[A-| ] generalized aerodynamic matrix 

{ dA } area, Eq. (6) 

{ F 0 } time Independent nonaerodynamlc 

force vector of a group of blades 

fj modal frequency of the blade In jth 

mode, HZ 

also blade Index 

f Kg J generalized stiffness matrix of a group 

of blades 

[K-j] generalized stiffness matrix of 1th 

blade In a group 


stiffness matrix of a group of 
blades, Eq. (1) 

centrifugal softening matrix of a group 
of blades 

mach number of freestream 

physical mass matrix 

generalized mass matrix of a group of 
blades 

generalized mass matrix of Jth blade 

generalized mass of the 1th blade In jth 
mode 

number of blade groups 

number of blades In a group 

number of modes of a group 

number of aerodynamic panels In a group 

number of modes of the 1th blade In a 
group 

total number of panel edge radii In a 
group 

stiffness matrix, Eq. (8); also loads 
on the blade group, Eq. (5) 


{PUu 0 })} 

steady state aerodynamic nodal 
force vector of a group of blades 


{q} 

generalized coordinate vector of 

a group 

1Q 0 > 

amplitude of { q } , Eq. (7) 


q U 

generalized coordinate jth mode of the 
1th blade In a group 

qi 

generalized coordinate vector of 
blade of a group 

1th 

{u 0 ) 

steady state deflection 


[W] 

matrix of normal velocities, Eq. 

0>) 

0Q.75R 

blade pitch angle at 3/4 radius 



[K( {u } ) ] 
[K s ] 

M 

[M] 

[M g 0 

[Mj] 

NG 

NB 

NMG 

NP 

NM1 

NRT 

[P] 
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y eigenvalue, Eq. (9) 

{ Au( t ) } vibratory deflection vector at grid 
points measured from steady state 


affect flutter and forced response; and (2) Inten 
tlonal mlstunlng may have potential as a passive 
control to Increase flutter speed and minimize 
system response. 


m 


v 

v 

a r 

°G1 

W U 


normal displacement, Eg. (6) 

damping radio In jth mode of 1th blade 
In a group 

Imaginary part of eigenvalue 

real part of eigenvalue 

Interblade phase angle between the 
blades of a group 

Intergroup phase angle between like 
blades from group to group, 

1 = 0,1,2,. ..(NG-1) 

model matrix of 1th blade In a group 

frequency of the Jth mode of the 1th 
blade 


Introduction 


Several questions arise concerning the role 
of mlstunlng In propfan flutter. Specifically: 

(1) Does the presence of sweep and strong bending 
torsion coupling between the blade structural modes 
Influence the role of mlstunlng? (2) Does mls- 
tunlng have enough potential to be used as a pas- 
sive control to alleviate flutter? (3) Does the 
presence of mlstunlng have any Influence on blade 
loads? To answer the first two questions par- 
tially, an analytical and experimental Investiga- 
tion was conducted. The analytical effort Involves 
an extension of the flutter model described In 
Ref. 5 to Include both structural and aerodynamic 
mlstunlng. The structural model for each blade Is 
defined by blade geometry, normal modes, frequen 
cles, and generalized masses. The aerodynamic 
model Is based on three-dimensional subsonic lift- 
ing surface theory. 8 The models account for 
differences In blade properties such as geometry, 
pitch angle, sweep, modes, frequencies, and damp- 
ing. The type of mlstunlng In the analytical 
flutter model Is arbitrary. 


Research on propfan aeroelastlcl ty Is a con- 
tinuing effort at NASA Lewis Research Center. 

Until now the analytical work has Involved propfans 
with blades of assumed Identical properties. This 
paper continues that work for propfans with blades 
of different properties. 

A flutter study of a propfan wind tunnel model 
with ten titanium blades under transonic tip condi- 
tions was reported In Ref. 1. Analytical compari- 
sons to the flutter data, presented In Ref. 1, were 
reported In Refs. 2 and 3. Additional flutter data 
from a model with eight composite material blades 
In subsonic flow was reported In Ref. 4. 


A wind tunnel experiment was performed. An 
elght-bladed mlstuned rotor was constructed by 
alternately mounting two sets of composite blades 
of two geometrically and materially similar rotors. 
The first set of four blades, SR3C-X2, are from 
the rotor described In Ref. 4. The second set of 
blades, SR3C-3, are Identical to the first set 
except for the ply angles of the composite mate 
rial. Because of the ply angle difference, the 
experimental rotor had a particular type of alter- 
nate blade mlstunlng which Is characterized by the 
differences In blade frequencies, mode shapes, 
steady state geometry, damping, and steady and 
unsteady air loads. 


A more refined analysis than those of Refs. 1 
to 3 Is described In Ref. 5. The analysis makes 
use of recent developments In modeling of composite 
material, 8 vibration analysis under rotating con- 
ditions Including geometric nonl Inearltles , ' and 
both steady and unsteady three-dimensional subsonic 
cascade aerodynamic models. 8 Additionally, 

Ref. 5 correlates predicted flutter results with 
the measured data reported In Ref. 4 for a propfan 
with eight Identical composite blades. 

In the analytical aeroelastlc models described 
In Refs. 1 to 3 and 5, It Is assumed that all the 
blades have Identical structural and aerodynamic 
properties. Consequently, there Is no coupling 
between the Interblade phase angle modes. Then, 
each mode can be analyzed Independently for flutter 
and forced response of the rotor. However, In 
practical propfans properties differ from blade to 
blade because of manufacturing limitations. The 
presence of small differences between the blades 
Is known as mlstunlng. Mlstunlng couples the 
Interb.ade phase angle modes, and Introduces com- 
plexity In the dynamic analysis. 

Understanding the effects of blade mlstunlng 
on vibration, flutter and forced response of turbo- 
machinery blade-disk assemblies has received con- 
siderable attention and Is a current research topic 
In turbomachinery aeroelastlclty . 9 ~1 4 The doml- 
nar 4 reasons fo r :nis attention are: (1) It does 


Flutter speed, flutter frequencies, and 
flutter Interblade phase angles were measured In 
the wind tunnel at various blade pitch angles, Mach 
numbers, and rotational speeds. Also measured 
during the test were blade strain amplitudes which 
will not be addressed In this paper. 

This paper presents the analytical formulation 
of the flutter model, describes the wind tunnel 
model and experiment, outlines the application of 
the formulation to the wind tunnel model, and com 
pares calculated and measured In-vacuum nonrotating 
frequencies and modes. Also compared In the paper 
are calculated and measured flutter speeds, flutter 
frequencies, and Interblade phase angles. Finally, 
the potential for alternate frequency mlstunlng as 
a passive flutter control Is assessed. 

Analytical Formulation 

The aeroelastlc analysis of propfans Is Inher 
ently nonlinear because of rotation, and requires a 
geometric nonlinear theory of elasticity. By using 
a finite element structural model and three 
dimensional steady and unsteady cascade aerodynamic 
models, an aeroelastlc model for a propfan with a 
rigid hub and Identical blades was developed In 
Ref. 5. That model has been generalized for a 
propfan with a rigid hub and with structurally and 
aerodynaml cal 1 y mlstuned blades. Since the basir 
details were given in Ref. S, only an outline oi 
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the development and of special modifications which 
are essential to allow mlstunlng are discussed 
here. 

Equations of Motion of a Hlstuned Propfan 

An aeroelastlc model with frequency mlstunlng 
was developed In Refs. 10 to 13 by exploiting the 
periodic structure of the aerodynamic model for a 
rotor with geometrically Identical blades by using 
two-dimensional unsteady cascade aerodynamic 
theory. Aerodynamic mlstunlng with supersonic two- 
dimensional unsteady cascade aerodynamic theory was 
addressed In Ref. 14. Since the present formula- 
tion does not distinguish between aerodynamic and 
structural mlstunlng and uses the normal modes of 
a rotating structure In conjunction with three- 
dimensional unsteady aerodynamic theory, a slightly 
different and more general approach than that used 
In Ref. 10 Is adopted. 


The next step In the formulation of the aero 
elastic model by the modal method Is to generalize 
the single blade equation of Ref. 5 to NB blades 
In a group. Assuming simple harmonic motion and 
following the procedure In Ref. 10, the aeroelastlc 
equations for NB arbitrarily mlstuned blades In a 
group are 

fH g J {q> ♦ fK g J {q} = [A ] ] {q} (3) 

Again, the form of Eq. (3) Is the same as the cor- 
responding one of Ref. 5, but the orders of the 
matrices and of the generalized coordinate vector 
are different. The new definitions are as follows: 

rt.,1 i 


In the present formulation, the propfan Is 
assumed to have NG Identical groups of blades 
symmetrically distributed about the disk. Each 
group contains NB blades (disjoint surfaces) 
which need have no special geometric or spatial 
relationship to each other. A schematic for an 
elght-bladed rotor with four blade groups Is Illus- 
trated In Fig. 1. A tuned rotor, then, corresponds 
to NB = 1, In which case NG Is the total number 
of blades. Conversely, If NG = 1 , then each of 
the blades, NB In number, Is treated Individually. 
In general. It Is assumed that the groups vibrate 
with an Identical motion, but with a constant 
Intergroup phase angle, which may take any of the 
values = 2ir1 /NG , 1 = 0,1 NG - 1. Gener- 

alizing the steady state single blade equation In 
Ref. 5, the corresponding equation for any group Is 

[[K s ] ♦ [K({U g ])]] <u o ) = {P({u o m ♦ { F o > (1) 
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The form of Eq. (1) Is the same as the correspond- 
ing one In Ref. 5, but the orders of the matrices 
and of the displacement vector are different. The 
stiffness and force matrices In Eq. (1) represent 
the entire blade group. For example, for a propfan 
with eight blades (four groups with two blades In 
each group) 228 grid points for each blade, and six 
degrees of freedom for each grid point, the total 
number of degrees of freedom of Eq. (1) Is 2736. 
Furthermore, the equation Is nonlinear Involving 
large deflections, centrifugal and steady state 
airloads. The solution requires substantial com- 
puter time. To reduce that time, the steady state 
aerodynamic loads have been neglected herein. (The 
effect of these steady airloads for a tuned rotor 
Is addressed In Ref. 5.) With this assumption, the 
steady aerodynamic coupling between the blades In a 
group Is eliminated from the flutter analysis. 
Consequently, the steady state deflections and the 
differential stiffness of each blade In a group can 
be calculated Independently by using Eq. (1) with- 
out steady airloads. The procedure for these cal- 
culations Is the same as that described In Ref. 5. 
The vibration modes and frequencies of each blade 
In the group are calculated Independently by using 
the following equation, which Is also discussed In 
Ref. 6. 
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The coefficients j , j , and Cll 
are the generalized mass, the corresponding fre- 
quency, and the corresponding damping ratio, 
respectively, of the Jth normal mode of the 1th 
blade In a group. The suffix NM1 represents the 
number of blade normal modes considered for the 1 th 
blade In the flutter analysis. The quantities 
[ 4 ^ ] and {q^} are the normal mode matrix and 
generalized coordinate vector for the ith blade In 
a group. The details for the development of the 
aerodynamic matrix [A^ ] are given In the follow- 
ing section. 

Unsteady Aerodynamic Model 

The generalized aerodynamic force matrix In 
Eq. ( 3 ) has been computed with the three- 
dimensional compressible lifting surface theory 
described In Refs. 8 , 15, and 16. Since the 
details of the scheme are given In those refer- 
ences, only an outline of the method and the spe- 
cial mlstunlng modifications will be discussed 
herein. 

The numerical method Is based on an assumed 
linearization of the fluid motion about a uniform 
steady flow. The disturbance generated by the 
rotor blades Is presumed to be simple harmonic In 
time In a reference frame rotating with the blades. 
Each wake Is assumed to be rigid and to lie In the 
helical surface swept out by the blade trailing 
edge. These assumptions result In a linear Inte- 
gral equation relating the normal velocity of the 
surface of the blades to the load (differential 
pressure) acting on the blades. The Integral equa- 
tion Is solved by a piecewise constant load panel 
method . 

For simplicity, the effects of hubs and 
nacelles are Ignored, as Is blade thickness. The 
blades are, in effect, treated as loaded helical 
surfaces with no thickness. As In any linear aero- 
dynamic theory, the unsteady loads Induced by blade 
vibration are mathematically decoupled from the 
steady loads arising from blade Incidence. 

The use of linear aerodynamic theory limits 
the applicability of the method to lightly loaded 
conditions (unstalled) and to Mach numbers at which 
transonic effects at the blade tips are not signi- 
ficant. In addition, the rigid wake approximation 
breaks down at low advance ratios, where Induced 
velocities may have a strong Influence on the wake 
structure (which In turn may significantly modify 
the blade loading.) Of course, one usually does 
not know aprlorl whether such effects will be 
present at a given operating point, so some caution 
In Interpreting the results Is necessary. 

In the original code. Ref. 5, on which the 
present method Is based, It Is assumed that the 
rotor contains NG equally spaced Identical blades, 
which vibrate with Identical motion but with a 
constant Interblade phase angle, one of the set 
2ir1/NG, 1 - 0,1, ...NG - 1. In the present version, 
mlstunlng Is accounted for by the simple expedient 
of redefining a blade as a group of NB Identical 
blades. Then these groups are symmetrically dis- 
tributed about the disk. 

The geometry and vibration mode data Is 
usually obtained by solving Eqs . (1) and (2) with 
NASI RAN. Each blade Is translated aiung ar^! 


rotated about the axis of rotation, and rotated 
about Its pitch change axis, so that It lies In the 
desired position and orientation. The group Is 
paneled by placing a fixed number, NXP, of quad- 
rilateral panels on each radial Interval 
rt(j) < r < rt( j * 1), J = 1,... NRT (unless 
rt(J + 1) < rt(j), which signals a Jump from the 
tip of one blade to the root of the next). Thus, 
except for the deletion of tip-root connections, 
the group Is effectively treated as a single sur- 
face with NP = NXP* (NRT - NB) panels. 

The definition of group vibration modes 
requires special mention. Essentially, one group 
mode consists of one blade In the group vibrating 
In one of Its natural modes while all other blades 
In the group remain fixed. Suppose the jth blade 
In a group Is assigned NMj vibration modes. The 
group as a whole, then has NMG * NM1 + NM2 + ... 

* NMNB modes. The group mode NM1 +1, for 
example, corresponds to the second blade In the 
group vibrating In Its first mode, while all 
remaining blades are stationary. 

In general, let W(n,k) denote the normal 
velocity at the nth control point In the kth group 
mode. This array Is Initialized to zero and then 
loaded sequentially with the NM1 values of first 
blade modes [<$>]], at control points on the 
first blade, the NM2 values of second blade 
modes, [$ 2 ], control points on the second 
blade, and so on. 

The aerodynamic loads can then be found by 
solving the linear algebra problem: 

[A] [P] = [W] (5) 

where [A] Is the (NP,NP) matrix of Influence coef- 
ficients, [W] Is the (NP, NMG) matrix of normal 
velocities, and [P] Is the (NP, NMG) array of loads 
on the blade group. The Influence coefficient 
matrix contains the Interference between blade 
groups and, therefore, depends on the Intergroup 
phase angle (as well as the vibratory frequency.) 
Thus, separate solutions are required for every 
Intergroup phase angle and frequency of Interest. 

Once the load distribution has been found, 
the generalized force matrix can be determined by 
simple numerical Integration over the blade group: 

NP 

[A, (n.m) ] = £ [ P( J ,m) ] [4(J,n)] {dA( j ) > ( 6 ) 

' J-l 

where P(j,m) Is the load on the Jth panel (with 
area dA(j)), associated with the mth group mode, 
and *(J,n) Is the normal displacement of the sur- 
face at the control point of the Jth panel In the 
nth mode. 

Solution 

For simple harmonic motion 

{q} = lV e1ut (7) 

Eq. (3) leads for the following: 

[P] (Qo) ■ irfHgJ ( 8 ) 
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where 


[P] * [KgJ - [A!] 

= 1w = v + 1v 


Flutter occurs when y = 0. The eigenvalue prob- 
lem, Eq. (8), Is solved for a given rotational 
speed and blade pitch angle and by using the proce- 
dure and the computer program, ASTR0P3, described 
In Ref . 5, for each of the possible Intergroup 
phase angles. Then, the flutter Intergroup phase 
angle Is the one for which the flutter Mach number 
Is the lowest. The Interblade phase angle In the 
group Is determined from the corresponding eigen- 
vector {q 0 }. From the calculated Intergroup 
and Interblade phase angles, the Interblade phase 
angles of the entire rotor are determined. 


Description of Experiment 


The experiment was conducted In the NASA Lewis 
8-by-6 ft. (2.44-by-l .83-m) wind tunnel. Test con- 
ditions Included tunnel Mach numbers from 0.36 to 
0.75 and rotor speeds up to 8000 rpm. The propfan 
models were mounted on a single rotation. Isolated 
nacelle test rig and driven by an air turbine. 

Blade mounted strain gauges provided blade vibra- 
tory strain signals. Eight blades were mounted In 
a hub which can be considered rigid. Fig. 2 shows 
the test rig Installed In the wind tunnel. 


Two existing propfan research models, SR3C-X2 
and SR3C-3, of 0.62 m (2 ft) diameter were used 
for this experiment. The blades were molded from 
graphlte-ply/epoxy-matrlx laminated material. Both 
models had 80 percent of the plys In the blade 
pitch axis (0°) direction, as shown In Fig. 3. The 
remaining plys were distributed at the ±22.5° 
directions for SR3C-X2, and at the ±45° directions 
for SR3C-3 . 


measured frequencies and mode shapes. These dif- 
ferences will be discussed later. 

RESULTS AND DISCUSSION 

The calculated and measured results for the 
SR3C-X2 and mixed rotors are presented In five 
parts: (1) measured nonrotating blade frequencies 

and mode shapes; (2) calculated blade rotating fre- 
quencies; (3) calculated root locus plots; (4) 
calculated damping plots; (5) calculated and meas- 
ured flutter results. Also, the calculated flutter 
characteristics for an alternately frequency mis- 
tuned SR3C-X2 rotor are presented. 

In all the calculations the steady state 
deflections, frequencies and mode shapes of the 
blades were obtained by using COSMIC NASTRAN. All 
the calculations are made by using the steady 
deflected position of the blade under centrifugal 
loads only. Only the first two modes for each 
blade were considered because, as will be shown 
later, the flutter mode Is primarily a combination 
of these two modes. Also for unsteady aerodynamic 
calculations each blade Is discretized with nine 
(radial) times eight (chordwlse) panels. 



The natural frequencies and mode shapes of 
the SR3C-X2 and -3 model differ because of the ply 
angle variation between the blades. Table I sum- 
marizes the measurements of the first two nonrota 
ting natural frequencies of the SR3C-X2 and -3 
blades for both rotor configurations. The first 
mode frequencies of both the blades are very close 
and were Insensitive to the variation of ply 
angles. However, the average second mode frequency 
of the SR3C-3 blade Is about 12 percent higher than 
that of the -X2 blade. 


The models were Individually wind tunnel 
tested prior to the present experiment. The 
SR3C-X2 model fluttered and the SR3C-3 did not 
within the operating conditions of the test (tunnel 
Mach = 0.9 and 10 000 rpm). The observed behavior 
of the two models differed because of frequency and 
mode shapes differences, which were, In turn, 
caused by the ply angle differences. Experimental 
results from the previous SR3C-X2 experiment were 
reported In Refs. 4 and 17. Analytical results 
were reported In Refs. 5 and 17 for the SR3C-X2 and 
In Ref. 17 for the SR3C-3. 

The present experiment was conducted In two 
parts, following the test procedure described In 
Ref. 4. First, the elght-bladed SR3C-X2 rotor was 
tested. This part Is a repetition of the test 
reported In Ref. 4 and was done because of changes 
In measured blade natural frequencies and mode 
shapes after the previous flutter test. The aver- 
age blade first and second mode frequencies (non- 
rotating) dropped by 5 and 1 percent, respectively, 
because of the blade material degradation from the 
large vibratory motion during flutter. Second, a 
mixed rotor, SR3C- X2/SR3C- 3 , was tested to obtain 
the flutter boundary of a configuration of blades 
with different frequencies and mode shapes. This 
rotor was constructed by alternately positioning 
four SR3C-X2 blades and four SR3C-3 blades. The 
SR3C-X2 rotor can be considered tuned and the mixed 
rotor mistuned based on A he differences between 


Figure 4 shows the nonrotating SR3C-X2 and -3 
blade displacement contours for the first two 
natural modes. Here, there Is a difference between 
the Inclination of the first mode contour lines of 
the blades. The greater Inclination of the SR3C-X2 
contour lines Indicates a greater torsional motion 
of the blade. The contour lines of the second 
modes of both blades are similar and show tip tor- 
sion. Thus, there are differences In the mode 
shapes of the SR3C-X2 and -3 blades, and the dif- 
ference Is more significant for the first mode. 

The. calculated and measured mode shapes are In fair 
agreement for both the blade models. 

A variation around the rotor of natural fre- 
quencies existed for both the SR3C-X2 and the mixed 
rotors but a mode shape variation of significance 
existed only for the mixed rotor. However, no 
attempt was made to quantify the mode shape differ- 
ences. The variation of the nonrotating measured 
first and second natural frequencies around the 
rotors are shown In Fig. 5. For the SR3C-X2 rotor, 
the largest natural frequency differences occurred 
between blades 2 and 3 (first mode, 6.5 percent; 
second mode, 5.5 percent). Whereas, for the mixed 
rotor, the largest difference occurred between 
blades 3 and 4 (first mode, 2.7 percent; second 
mode, 14.3 percent). Hence, the frequency mistun 
Ing was higher In the first mode for the SR3C- X2 
rotor and higher In the second mode for the mixed 
rotor . 
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Blade Frequency Variation with Rotational Speed 

Figures 6(a) and (b) show the variation with 
rotational speed of the blade frequencies In vacuum 
for the SR3C-X2 blade and the SR3C-3 blade respec- 
tively. The nonrotating average measured frequen- 
cies from Table I are also shown. The results show 
an Increase In the first mode frequency with rota- 
tional speed, whereas, the second mode frequency 
Is relatively constant until 6000 rpm and then 
Increases slightly. This Is expected because the 
first mode Is primarily bending, which Is affected 
by rotational speed, whereas, the torsion mode Is 
relatively Insensitive to rotation. The results 
shown are for the 61 .2° blade pitch angle. 

Results, although not shown, for a 68.0° blade 
angle are very close and follow the same trends as 
In Figs. 6(a) and (b). This Indicates Insensitiv- 
ity of the frequencies to a blade angle variation 
of this order. It should be mentioned that there 
Is a problem In estimating the equivalent aniso- 
tropic material properties, because of the unknown 
blade material properties after fabrication. When 
the material manufacturers' specified properties 
of the composite material were used In calculating 
blade nonrotating frequencies, the frequencies were 
higher than those measured. Therefore, to obtain 
better agreement the material constants were 
adjusted . 

Calculated Root Locus Plots 

To Illustrate the effect of mlstunlng (which 
Is partly aerodynamic and structural) the calcula- 
ted real and Imaginary parts of the eigenvalues of 
the SR3C-X2 , SR3C-3 and mixed rotors were compared. 
The mlstunlng is due to the differences In blade 
steady state geometry, frequencies, and mode 
shapes. The steady state geometry of the blades 
In a group differs because the steady deflections 
due to centrifugal loads are different. The eigen- 
values are for all Interblade phase angles of the 
mode with lowest damping at the critical Hach num- 
ber and critical reduced frequency of the SR3C-X? 
rotor. The calculations were performed by treating 
the SR3C-X2 and SR3C-3 rotors as tuned and the 
other one as mlstuned with both structural and 
aerodynamic mlstunlng. Furthermore, the mixed 
rotor Is modeled as an Idealized alternately mls- 
tuned rotor -- four Identical blade pairs with two 
different blades In each pair. The effective damp- 
ing ratio of each Interblade phase angle mode Is 
given by the ratio of the real and Imaginary parts 
of that mode. 

The comparison of the root loci for the 
SR3C-X2 and SR3C-3 rotors (Fig. 7) reveals two 
Interesting points. First, the area of the approx- 
imate ellipse for the SR3C-X2 Is greater than that 
of the SR3C-3 , Indicating a stronger aerodynamic 
coupling between the blades of the SR3C-X2 rotor. 
Second, at the chosen Mach number the SR3C-X2 rotor 
Is neutrally stable and the SR3C-3 rotor Is stable. 
This difference In stability may be attributed to 
the differences In blade stiffness and mode shapes 
and to the degree of aerodynamic coupling between 
the normal modes of the blades. The blade stiff- 
ness and mode shapes are different because of the 
different ply angles of the blades. Thus, one can 
Infer that a laminated composite propfan blade can 
be tailored to maximize Its flutter speed by 
selecting the proper ply angles. A comparison of 
the eigenvalues (Fig. 7) for the SR3C-X2, SR3C-3, 


and mixed rotors shows that mixing the blades sig- 
nificantly affected the eigenvalues and resulted 
In a rotor with a greater damping than the lowest 
damped mode of either pure rotor. 

The aerodynamic coupling between the modes of 
neighboring blades Is virtually eliminated In the 
mixed rotor. The effective damping of some Inter- 
blade phase angle modes Is Increased while that of 
others Is decreased. Additionally, the root locus 
plot of the mixed rotor Is divided Into high and 
low frequency compact groups such that the effec- 
tive damping of each Interblade phase angle of each 
group Is almost the same. The high frequency group 
of eigenvalues lies Inside the SR3C-X2 ellipse and 
the other group Inside the SR3C-3 ellipse. It Is 
the effect of the aerodynamic coupling of the modes 
of both the SR3C-X2 and -3 blades In the mixed 
rotor that give the low and high frequency eigen- 
value clusters. The modes of the high frequency 
cluster Involve motion primarily of the SR3C-X2 
blades and the modes of the low frequency cluster 
Involve motion primarily of the -3 blades. The 
close spacing of the eigenvalues In each group 
shows that the aerodynamic coupling between neigh- 
boring blades (the cascade effect) has been nulli- 
fied. Note that both the 180 and 225° Interblade 
phase angle modes are neutrally stable for the 
SR3C-X2 rotor. But, the flutter frequencies are 
different. 

To see the effect of Mach number on root loci, 
the calculated root locus plots (not shown) for the 
three rotors were also compared at the critical 
Mach number and reduced frequency of the mixed 
rotor. For this case the root locus of the SR3C-X2 
rotor shown In Fig. 7 moved to the right, such that 
four of the Interblade phase angle modes are In the 
unstable region. But, the root locus of the SR3C-3 
rotor moved very little. Also, the cluster of four 
higher frequency eigenvalues of the mixed rotor 
remained centered In the SR3C-X2 ellipse and, 
hence, moved close to the neutral stability line. 
But, the other cluster, remaining roughly centered 
In the SR3C-3 locus, moved very little. This shows 
that the change In damping of the modes In the high 
frequency cluster Is more sensitive than the low 
frequency cluster to the Increase In Mach number. 
The effect of Mach number on damping will be 
further discussed In Fig. 8. 

In summary, the mlstunlng has Increased the 
damping of the lowest damped modes of both the pure 
rotors and has nullified the cascade effects. 

These are the same trends that have been found for 
turbofan rotor In Ref. 10, although these propfans 
have fewer blades, lower gap-to-chord and mass 
ratios, and mode shape and aerodynamic mlstuning. 

Calculated Damping 

Figure 8 shows the variation of the real part 
of the eigenvalue (which Is propor t Iona 1 to damp 
Ing) with freestream Mach number at Q - 6320 rpm 
and Bo. 75R = 61.2° for the lowest damped mode of 
the SR3C-X2, SR3C-3, and the mixed rotors. Damping 
curves for the 180 and 226° interblade phase angle 
modes of the SR3C-X2 rotor are shown because both 
the modes become unstable at the same Mach number. 
Note that there are two curves for the mixed 
rotor: one for the low frequency cluster and the 

other for the high frequency cluster of Mg. 7. 
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If the mixed rotor Is considered to be 
developed by perturbation from the SR3C-X2 rotor, 
the damping of the least stable mode has Increased 
for all the Mach numbers shown. If the mixed 
rotor Is constdered to be developed by 
perturbation from the SR3C-3 rotor the change In 
damping depends on the Mach number. For Mach 
numbers between 0.40 and 0.55 the damping has 
Increased. Beyond Mach numbers of 0.55 the 
damping of the high frequency mode of the mixed 
rotor Is greater than that of the low frequency 
mode and Is less than that of the SR3C-3 rotor. 

The change In the damping from each of the pure 
rotors Is Illustrated by the arrows In the 
figure. Thus, the mixing of the two sets of 
blades has a beneficial effect on damping with 
respect to both of the pure rotors between Mach 
numbers of 0.40 to 0.55. Whereas, above 0.55 Mach 
number the mixing has a beneficial effect on the 
damping with respect to the SR3C-X2 rotor but an 
adverse effect on the damping with respect to the 
SR3C-3 rotor. 

Flutter Results 

Experimental and calculated flutter 
boundaries are correlated In Fig. 9 for 
Bq . 75R = 61.2°. Results for the SR3C-X2 and 
mixed rotors are Included. The calculated 
boundary for each rotor was obtained with the 
calculated modes and frequencies, except that the 
measured second mode frequency was substituted for 
the calculated one. The calculated flutter Mach 
numbers for the SR3C-X2 are less than the measured 
ones for all rotational speeds. The agreement 
between theory and experiment Is expected to be 
better If the effects of steady state aerodynamics 
and the structural damping are Included In the 
analysis (see Ref. 5 for detailed discussion). 

The agreement for the mixed rotor Is better (with 
the present theory), but would become 
unconservative If steady-state aerodynamic loads 
and structural damping were Included In the theory. 

The correlation of measured and calculated 
flutter frequencies at a blade angle of 61.2° Is 
given In Figs. 10(a) and (b) for the SR3C-X2 and 
the mixed rotors, respectively. The flutter 
frequencies are shown on a Campbell diagram which 
Includes the calculated frequencies for the first 
two modes of the blades, as well as the 2P, 3P, 
and 4P harmonic order excitation lines of 
rotational speed. For both figures, the 
calculated flutter frequencies are higher than the 
measured frequencies by about 30 Hz, but the 
slopes of the analytical and measured flutter 
frequency lines are about the same. The fact that 
both the calculated and measured frequencies lie 
between the first and second normal mode 
frequencies shows that the flutter mode Involves a 
coupling between these two modes. 

Measured and calculated phase angles for the 
mixed rotor are shown In Table II for three 
conditions. The measured values were obtained 
from cross-spectra of blade strain gauge signals. 
The coherence values associated with the phase 
angles were not always good and only phase angles 
with coherence values greater than 90 percent are 
shown. The data shows that several Interblade and 
Intergroup phase angles are present during flutter 
and that there are no predominant Intergroup or 
Inters ade phase angles. The analysis also showed 
that several Interblade and Intergroup phase angle 


modes were very lightly damped, however, only the 
lowest damped mode Is shown In Table II. So quali- 
tatively the theory and experiment are In agree- 
ment, however, the lack of quantitative agreement 
may be due to the fact that the analytical model 
assumes Identical properties for all the blade 
groups while the actual rotor has differences 
between blade groups. For example, the differences 
of blade frequencies Is Illustrated In Fig. 4. 

The results shown In Figs. 9 and 10 and 
Table II were also calculated and measured for 
Bq, 75R ■ &8°. Although not shown, these results 
reveal similar trends and correlation between 
theory and experiment. However, the flutter speeds 
are greater by approximately 0.04 Mach number for 
the blade angle of 68°. In summary, the overall 
agreement between theory and experiment Is reason- 
ably good. 

Calculated Flutter Characteristics of an 
Alternately Frequency Mlstuned Rotor 

The flutter results presented In Fig. 7 for 
the SR3C-X2 rotor shows that the aerodynamic cas- 
cade effects Is significantly destabilizing. 
Quantitatively, the destabilizing effect Is approx- 
imately proportional to the length of the semi- 
major axis of an ellipse formed by Joining all the 
eight points. This destabilizing cascade effect 
was shown to be present, even for a four-bladed 
rotor In Ref. 5. Since the potential beneficial 
effect on damping of frequency mlstunlng In a 
turbofan rotor (Ref. 10) Is directly proportional 
to the adverse effect of aerodynamic cascade 
effects. It Is natural to ask whether a similar 
relationship between cascade and mlstunlng effects 
exists for propfan rotor. To answer this question, 
six configuration cases, as shown In Table III, of 
the SR3C-X2 rotor are considered. The calculated 
flutter Mach number, flutter frequencies and 
flutter Interblade phase angles are Included In 
this table. The rotor In Case 1 Is the reference 
rotor. In order to confirm the Importance of cas- 
cade aerodynamic effects, the results of eight, 
four, and one-bladed tuned rotors (Cases 1, 2, 
and 3) are compared. Evidently, the flutter Mach 
number decreases with Increasing number of blades. 
The decrease In flutter Mach number due to cascade 
effects with eight blades Is 0.10. 

In order to evaluate the sensitivity of 
flutter Mach number to blade frequency changes, all 
blade frequencies are decreased by 10 percent for 
the rotor In Case 4 and are Increased the same per- 
cent for the rotor In Case 5. The flutter Mach 
number for the Case 4 Is decreased from 0.53 to 
0.46 and for the Case 5 Is Increased from 0.53 to 
0.59, showing that the flutter Mach number Is sen- 
sitive to change In the blade frequencies. 

To establish the effects of alternate fre- 
quency mlstunlng on flutter, the frequencies of the 
even blades were Increased by 10 percent and those 
of the odd blades were decreased by the same per- 
cent (Case 6). Alternate blade frequency mlstunlng 
Is considered because the beneficial effects on 
flutter are "robust" for this kind of mlstunlng In 
turbofan rotors, see Refs. 10, 12, and 18. From 
this established relationship between the adverse 
cascade effects and beneficial alternate frequency 
mlstunlng effects (Ref. 10), the expected range of 
the flutter Mach number for the Case 6 Is between 
0.46 (Case 4) and 0.63 (Case 3). To be beneficial 
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as a passive flutter control, the flutter Mach 
number for Case 6 should be greater than 0.59 
(Case 5). But the calculated flutter Mach number 
for Case 6 Is 0.53. Hence, It Is concluded that 
alternate frequency mlstunlng does not have enough 
potential to be used as a passive control to 
Increase flutter speed of this kind of propfan. 
However, mlstunlng should be Included In the ana- 
lysis to assess the effective damping of each 
Interblade phase angle mode because the damping Is 
affected significantly by the presence of mlstunlng 
as seen In Fig. 7. The fact that the flutter Mach 
number for Case 6 Is lower than expected from 
turbofan studies may be due to the strong coupling 
between bending and torsion and blade sweep. 
Unfortunately, no attempt was made In this paper 
to Isolate the aerodynamic and structural effects 
of blade sweep even though the analytical model has 
the provision to account for It. Comparing the 
effects of blade alternate frequency mlstunlng on 
turbofan and propfan flutter. It appears that the 
beneficial effects of this kind of mlstunlng In 
propfan rotors are nullified by the adverse effects 
of the lower alternate blade frequencies. An 
Interesting point to be noted from Table III Is 
that the critical Interblade phase angle varies 
with the blade frequencies and number of blades. 

It also varies with blade pitch angle, as was shown 
In Ref. 5. 

CONCLUSIONS 

An analytical and experimental Investigation 
on the effect of both structural and aerodynamic 
mlstunlng on propfan subsonic flutter was con- 
ducted. An aeroelastlc model for the analysis of 
propfan flutter with mlstunlng was developed and 
was validated for selected cases by comparing the 
calculated and measured flutter characteristics of 
a propfan wind tunnel model with composite blades. 
It was found that the correlation between theoret- 
ical and experimental flutter Mach numbers Is good. 
But the calculated frequencies are higher than the 
measured by approximately about 3 to 6 percent. 

Both theory and experiment also showed that com- 
bined mode shape, frequency and aerodynamic mls- 
tunlng can have either a beneficial or adverse 
effect on blade damping depending on the Mach num- 
ber. Additional parametric results showed that 
alternate blade frequency mlstunlng does not have 
enough potential for It to be used as a passive 
control to alleviate flutter In propfans similar 
to the one studied herein. The beneficial effects 
of blade alternate frequency mlstunlng on flutter 
were nullified by the adverse effects of blade fre- 
quency drop due to alternate frequency mlstunlng. 
Finally, It Is Inferred from both theoretical and 
experimental results that a laminated composite 
propfan blade can be tailored to optimize Its 
flutter speed by selecting the proper ply angles. 
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TABLE I. - N0NR01 AT ING BLADE MEASURED FREQUENCIES 


Rotor 

configuration 

Blade 

type 

Number of 
blades 

Mode 1 

Mode 2 

Average 

Range 

Average 

Range 

SR3C-X2 

- X2 

8 

191 

185-197 

371 

362-383 

SR3C-X2/ 

- X2 

4 

189 

185-193 

369 

363-378 

SR3C-3 

- 3 

4 

188 

185-190 

415 

415-416 


TABLE II. - MEASURED AND CALCULATED PHASE ANGLES FOR THE 
SR3C- X2/SR3C- 3 ROTOR AT A BLADE ANGLE OF 61.2° 



Measured 3 

Calculated 



°r* 

deg 



a G> 

deg 


a r » 

0 G» 










deg 

deg 

Blades -> 

1-2 

3-4 

5-6 

7-8 

1-3 

3-5 

5-7 

7-1 



RPM 











5230 

154 

326 

121 

192 

324 

250 



_ 

155 

180 

6750 

138 

175 

— 

236 

211 

148 

— 

224 

155 

90 

7450 

139 

188 

— 

229 

205 

134 

— 

202 

146 

90 


a The phase angle of blades A-B Is that of B with respect to A. 


TABLE III. - FREQUENCY MISTUNING STUDY OF THE SR3C-X2 PROPFAN; BLADE ANGLE, 
61.2°, ROTATIONAL SPEED, 6320 RPM 


Case 

Number of 
blades 

Blade 

A frequency, 
percent 

Flutter conditions 

Freestream, 

Mach 

Frequency, 

Hz 

Phase 

Angle, <? r , 
deg 

1 . 

Tuned 

8 

0 

0.53 

294 

225 






302 

180 

2. 

Tuned 

4 

0 

0.60 

303 

180 

3. 

Tuned 

1 

0 

0.63 

306 

0 

4. 

Tuned 

8 

-10 

0.46 

277 

180 

5. 

Tuned 

8 

+10 

0.59 

317 

225 

6. 

Alternately 

8 

+10 

0.53 

284 

a 275 


mlstuned 








a The corresponding group phase angle Is 90°. 
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FIGURE 1.- BLADE GROUP SCHEMATIC FOR AN EIGHT-BLADED ROTOR. 
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FIGURE 2.- SR3C-X2 PROPFAN MODEL WIND TUNNEL INSTALLATION. 
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FIGURE 3.- BLADE PLY DIRECTIONS. 
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NON-ROTATING BLADE MEASURED AND CALCULATED MODE SHAPES. 





190 Hz (MODE 1) 
368 Hz (MODE 2) 



189 

365 

(A) SR3C-X2 ROTOR 


190 Hz (MODE 1) 
368 Hz (MODE 2) 



189 

365 


(B) SR3C-X2/SR3C-3 
ROTOR 

FIGURE 5.- NON-ROTATING BLADE FREQUENCY DISTRIBUTION 
AROUND THE ROTORS. 




0 2000 4000 6000 8000 

ROTATIONAL SPEED, rpm 

(B) SR3C-3 MODEL BLADE. 

FIGURE 6. - CALCULATED VARIATION OF BLADE NATURAL 
FREQUENCIES WITH ROTATION AT 61.2° BLADE ANGLE. 
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IMAGINARY PART OF EIGENVALUE 



FIGURE 7. - CALCULATED ROOT LOCUS PLOT OF THE LOWEST 
DAMPED MODE FOR THE SR3C-X2, SR3C-3. AND SR3C-X2/ 
SR3C-3 ROTORS AT 61.2° BLADE ANGLE, 6320 RPM AND 
0.528 FREESTREAM MACH NUMBER. 
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ROTOR 


O SR3C-X2 



FIGURE 8. - VARIATION OF CALCULATED DAMPING WITH MACH 
NUMBER FOR THE SR3C-X2 AND SR3C-X2/SR3C-3 ROTORS AT 
61.2° BLADE ANGLE AND 6320 RPM. 


17 



SR3C-X2 

ROTOR 


O 


z 

Cl 

01 


Q_ 

CO 


< 

o 

< 

5 


8000 

7500 

7000 

6500 

6000 

5500 

5000 

4500 


.3 




STABLE 


FIGURE 9. - ME, 
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FIGURE 10. - MEASURED AND CALCULATED FLUTTER FREQUENCIES 
AT 61.2 DEG BLADE ANGLE. 
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